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> ■ Abstract 
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, We discuss the possibility that in finite density QCD an anisotropic 

| phase is realized. This case might arise for quarks with different chem- 

ical potential and/or different masses. In this phase crystalline struc- 
tures may be formed. We consider this possibility and we describe, in 
f-^. | the context of an effective lagrangian, the corresponding phonons as 

^^O ■ the Nambu-Goldstone bosons associated to the breaking of the space 

. symmetries. 

m 

1 Introduction 

^p Hi In the last few years color superconductivity has received a lot of attention [1, 

2] . The possibility of formation of diquark condensates was pointed out since 
^ . the first days of QCD [3]. Naively one would expects that, due to asymptotic 

j> \ freedom, quarks at very high density would form a Fermi sphere of almost free 

fermions. However, Bardeen, Cooper and Schrieffer [4] proved that the Fermi 
surface of free fermions is unstable in presence of an attractive, arbitrary 
small, interaction. Since in QCD the gluon exchange in the 3 channel is 
attractive one expects the formation of a coherent state of particle/hole pairs 
(Cooper pairs). An easy way to understand the origin of this instability is 
to remember that, for free fermions, the Fermi energy distribution at zero 
temperature is given by f(E) = 9(fi — E) and therefore the maximum value 
of the energy (Fermi energy) is Ep = \i. Then consider the corresponding 
grand-potential, Q = E — /iiV, with \i = Ep- Adding or subtracting a 
particle (or adding a hole) to the Fermi surface does not change Q, since 

1 On leave from the Department of Physics of the University of Florence. 
2 Invited talk at "QCD@Work 2003 - International Workshop on QCD, Conversano, 
Italy, 14-18 June 2003". 



1 



Q — > (E ± — //(TV =b 1) = Q. We see that the Fermi sphere of free 
fermions is highly degenerate. This is the origin of the instability, because if 
we compare the grand-potential for adding two free particles or two particles 
bounded with a binding energy Eb, we find that the difference is given by 
il B — Q = —E B < 0. Since a bound state at the Fermi surface can be 
formed by an arbitrary small attractive interaction [5], it is energetically 
more favorable for fermions to pair and form condensates. 

From the previous considerations it is easy to understand why to realize 
superconductivity in ordinary matter is a difficult job. In fact, one needs an 
attractive interaction to overcome the repulsive Coulomb interaction among 
electrons, as for instance the one originating from phonon exchange for elec- 
trons in metals. On the other hand the interaction among two quarks in the 
channel 3 is attractive , making color superconductivity a very robust phe- 
nomenon. Notice also that, once taken into account the condensation effects, 
in the limit of very high density one can use asymptotic freedom to get exact 
results. For instance, it is possible to get an analytical expression for the gap 
[6]. In the asymptotic regime it is also possible to understand the structure 
of the condensate 3 . In fact, consider the matrix element 

(oiC^lo) (i) 

where a, (3 = 1,2,3 are color indices, a,b — 1,2 are spin indices and i,j = 
1, • • • ,N are flavor indices. Its color, spin and flavor structure is completely 
fixed by the following considerations: 

• antisymmetry in color indices (a, (3) in order to have attraction; 

• antisymmetry in spin indices (a, b) in order to get a spin zero conden- 
sate. The isotropic structure of the condensate is favored since it allows 
a better use of the Fermi surface; 

• given the structure in color and spin, Pauli principles requires antisym- 
metry in flavor indices. 

Since the momenta in a Cooper pair are opposite, as the spins of the quarks 
(the condensate has spin 0), it follows that the left(right)-handed quarks can 
pair only with left(right)-handed quarks. In the case of 3 flavors the favored 

3 At non asymptotic densities other possible condensates might be formed, see [7] 
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condensate is 



(0\€l^l\0) = -<0|«>> = A J2 (2) 

c=i 

This gives rise to the so-called color-flavor-locked (CFL) phase [8, 9]. The 
reason for the name is that simultaneous transformations in color and in 
flavor leave the condensate invariant. In fact, the symmetry breaking pattern 
turns out to be 

SU(3) C ® SU(3) L <g> SU(3) R <g> U(1) B -> SU(3) C+L+R 

where SU(3) C+L+R is the diagonal subgroup of the three SU(3) groups. This 
is the typical situation when the chemical potential is much bigger than the 
quark masses m u , and m s (here the masses to be considered are in prin- 
ciple density depending) . However we may ask what happens decreasing the 
chemical potential. At intermediate densities we have no more the support 
of asymptotic freedom, but all the model calculations show that one still has 
a sizeable color condensation. In particular if the chemical potential /i is 
much less than the strange quark mass one expects that the strange quark 
decouples, and the corresponding condensate should be 

<0|Cl^ l |0> = Ae Q/ % (3) 

In fact, due to the antisymmetry in color the condensate must necessarily 
choose a direction in color space. Notice that now the symmetry breaking 
pattern is completely different from the three-flavor case. In fact, we have 

SU(3) C <g> SU(2) L <g> SU(2) R ® U(1) B -> 

-> SU(2) C <g> SU(2) L <g> SU(2) R <g> U(1) B 

It is natural to ask what happens in the intermediate region of fi. It turns 
out that the interesting case is for /i « M s 2 /A. To understand this point let 
us consider the case of two fermions: one massive, mi = M s and the other 
one massless, rri2 = 0, at the same chemical potential /i. The Fermi momenta 
are of course different 

p Fl = vV - M*, p F2 = /I (4) 
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The grand potential for the two unpaired fermions is (the factor 2 comes 
from the spin degree of freedom) 



f PF 2 d 3 v 

In order to pair the two fermions must reach some common momentum p^ mm , 
and the corresponding grand potential is 



rPcomm 

l>a " ' = 2 Jo : " rr ^ 



d 3 p 



(2k) 



"pair. = 2 / -^-[ Jj? + M*-» 



where the last term is the energy necessary for the condensation of a fermion 
pair [10]. The common momentum Pc 0mm can be determined by minimizing 
r2p a i r . with respect to Pc 0mm , with the result 

F M s 
Pcomm = » ~ (7) 

It is now easy to evaluate the difference n unpair . — f2 pair . at the order , with 
the result 

^unpair. " ^pair. « ^ (M* - 4A 2 /i 2 ) (8) 

We see that in order to have condensation the condition 

M 2 

»>2K < 9 > 

must be realized. The problem of one massless and one massive flavor has 
been studied in ref. [11]. However, one can simulate this situation by taking 
two massless quarks with different chemical potentials, which is equivalent 
to have two different Fermi spheres. The big advantage here is that one can 
use a study made by Larkin and Ovchinnikov [12] and Fulde and Ferrel [13]. 
These authors studied the case of a ferromagnetic alloy with paramagnetic 
impurities. The impurities produce a magnetic field which, acting upon the 
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electron spins, gives rise to a different chemical potential for the two popu- 
lations of electrons. It turns out that it might be energetically favorable to 
pair fermions which are close to their respective Fermi surface (LOFF phase). 
However, since the Fermi momenta are different, the Cooper pair cannot have 
zero momentum and there is a breaking of translational and rotational invari- 
ance. Therefore, a crystalline phase can be formed. The previous situation 
is very difficult to be realized experimentally, but there have been claims 
of observation of this phase in heavy-fermion superconductors [14] and in 
quasi-two dimensional layered organic superconductors [15]. The authors of 
ref. [16] have extended the calculation of ref. [12, 13] to the case of two-flavor 
QCD and we will review here their results. Also, since the LOFF phase can 
give rise to crystalline structures, phonons are expected. We will also discuss 
the effective lagrangians for the phonons in different crystalline phases and 
show how to evaluate the parameters characterizing them, as the velocity of 
propagation. Finally we will consider an astrophysical application. 



According to the authors of refs. [12, 13] when fermions belong to two differ- 
ent Fermi spheres, they may prefer to pair staying as much as possible close 
to their own Fermi surface. When they are sitting exactly at the surface, 
the pairing is as shown in Fig. 1. We see that the total momentum of the 



Figure 1: Pairing of fermions belonging to two Fermi spheres of different 
radii according to LOFF. 

pair is pi + p2 = 2g and, as we shall show, \q\ is fixed variationally whereas 



2 The LOFF phase 
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the direction of q is chosen spontaneously. Since the total momentum of the 
pair is not zero the condensate breaks rotational and translational invariance. 
The simplest form of the condensate compatible with this breaking is just a 
simple plane wave (more complicated functions will be considered later) 

Oj{x)jP(x)) « Ae 2! ^ (10) 

It should also be noticed that the pairs use much less of the Fermi surface than 
they do in the BCS case. In fact, in the case considered in Fig. 1 the fermions 
can pair only if they belong to the circles in the figure. More generally there 
is a quite large region in momentum space (the so called blocking region) 
which is excluded from the pairing. This leads to a condensate smaller than 
the BCS one. 

Before discussing the LOFF case let us review the gap equation for the 
BCS condensate. We have said that the condensation phenomenon is the 
key feature of a degenerate Fermi gas with attractive interactions. Once one 
takes into account the condensation the physics can be described using the 
Landau's idea of quasi-particles. In this context quasi-particles are nothing 
but fermionic excitations around the Fermi surface described by the following 
dispersion relation 

e(p,A BCS ) = yJ? + A BCS (11) 

with 

Z = E(p)-»a?^ .(p-p F )=v F -(p-p F ) (12) 

Op p=p F 

The quantities vp and {p — Pf) are called the Fermi velocity and the resid- 
ual momentum respectively. A easy way to understand how the concept of 
quasi-particles comes about in this context is to study the gap equation at 
finite temperature. For simplicity let us consider the case of a four-fermi 
interaction. The euclidean gap equation is given by 

/ d'p 1 

9 J (27r)V! + b1 2 + A| cs 1 6) 

From this expression it is easy to get the gap equation at finite temperature. 
We need only to convert the integral over p 4 into a sum over the Matsubara 
frequencies 



J3 +°° 1 

i dp 1 

= // ' J^f „^ ((2n + PF + e*(p, A BCS ) { " ] 



n=—oo 
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Performing the sum we get 



, 9 f d 3 p l-n u -n d 

2 J (2tt)3 e(p,A BCS ) 1 } 

Here n u and are the finite-temperature distribution functions for the ex- 
citations (quasi-particles) corresponding to the original pairing fermions 

n ^ = n d = e ^ ABC ] )/T - i (16) 

At zero temperature (n u = rid — ► 0) we find (restricting the integration to a 
shell around the Fermi surface) 

vldi 1 



# r dttpjr 
2 J (2tt 



r ) 3 + a| C5 

In the limit of weak coupling we get 

A BCS ~2Ze- 2 /^ (18) 

where £ is a cutoff and 

P=# (19) 

is the density of states at the Fermi surface. This shows that decreasing 
the density of the states the condensate decreases exponentially. From a 
phenomenological point of view, one determines the coupling g requiring 
that the same four-fermi interaction at zero temperature and density, gives 
rise to a constituent mass of the order of 400 MeV. From this requirement, 
using values for ii w 400 -j- 500 MeV (interesting for the physics of compact 
stellar objects), one obtains values of A in the range 20 -j- 100 MeV . Since 
at very high density it is possible to use perturbative QCD one can evaluate 
the gap from first principles [17]. The result is 

A w bixgle-^l^ 9 * (20) 

It is interesting to notice that from Nambu-Jona Lasinio type of models one 
would expect a behavior of the type exp(— c/g 2 ) rather than exp(— c/g). This 
is due to an extra infrared singularity from the gluon propagator. Although 
this result is strictly valid only at extremely high densities, if extrapolated 
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down to densities corresponding to /i ~ 400 -r- 500 MeV, one finds again 
A « 20 4- 100 MeF. 

Let us now consider the LOFF case. For two fermions at different densities 
we have an extra term in the hamiltonian which can be written as 

Hj = -5fxa 3 (21) 

where, in the original LOFF papers [12, 13] 5/i is proportional to the magnetic 
field due to the impurities, whereas in the actual case 5[i = (/ii — /U 2 )/2 and 
cr 3 is a Pauli matrix acting on the two fermion space. According to refs. 
[12, 13] this favors the formation of pairs with momenta 

p 1 = k + q, p 2 = —k + q (22) 

We will discuss in detail the case of a single plane wave (see eq. (10)). The 
interaction term of eq. (21) gives rise to a shift in £ (see eq. (12)) due both 
to the non-zero momentum of the pair and to the different chemical potential 

£ = E(p) - n -> E(±k + q)- f i T 5fi^£Tfi (23) 

with 

ji = 5/i — vf • q (24) 

Here we have assumed S/i /i (with \i = (/xi + 7/2) /2) allowing us to expand 
E at the first order in q (see Fig. 1). The gap equation has the same formal 
expression as in eq. (15) for the BCS case 



1/ 



d 3 p l-n u -n d 

(25) 



(2vr)3 e(p,A) 



but now n u 7^ n d 



n ^ d e (e(p,A)±H)/T + 1 ( 26 ) 
where A is the LOFF gap. In the limit of zero temperature we obtain 

^f/lf^A)' 1 -^-^-^^" < 27 > 

The two step functions can be interpreted saying that at zero temperature 
there is no pairing when e(p, A) < This inequality defines the so called 
blocking region. The effect is to inhibit part of the Fermi surface to the 
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pairing giving rise a to a smaller condensate with respect to the BCS case 
where all the surface is used. 

We are now in the position to show that increasing 5 /a from zero we have 
first the BCS phase. Then at 5fi = 8li\ there is a first order transition to the 
LOFF phase [12, 16], and at 6 fx = 5fi2 > there is a second order phase 
transition to the normal phase (with zero gap) [12, 16]. We start comparing 
the grand potential in the BCS phase to the one in the normal phase. Their 
difference is given by 

2 

^BCS — ^normal = — T~9 i^BCS ~ 2<5/i 2 ) (28) 

47T Vp 

where the first term comes from the energy necessary to the BCS conden- 
sation (compare with eq. (6)), whereas the last term arises from the grand 
potential of two free fermions with different chemical potential. We recall 
also that for massless fermions pp = n and Vp — 1. We have again assumed 
5fi <C li. This implies that there should be a first order phase transition from 
the BCS to the normal phase at 811 = Abcs/V% [18], since the BCS gap 
does not depend on Sli. The situation is represented in Fig. 2. In order to 




Figure 2: The grand potential and the condensate of the BCS and LOFF 
phases vs. Sli. 

compare with the LOFF phase we will now expand the gap equation around 
the point A = (Ginzburg-Landau expansion) in order to explore the possi- 
bility of a second order phase transition. Using the gap equation for the BCS 
phase in the first term on the right-hand side of eq. (27) and integrating the 
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other two terms in £ we get 



9 p% log ^cs = jA f^^m (29) 



2tt 2 v f to A 2tt 2 v f J 4tt A 

where 

C(9) = ^{5p-qv F cos#) 2 - A 2 (30) 

For A — > we get easily 



log 



2<5/i 



i/ log(l-H), , = ^ (31) 



This expression is valid for 5p smaller than the value Sp 2 at which A = 0. 
Therefore the right-hand side must reach a minimum at 5p = 5p 2 . The 
minimum is fixed by the condition 

- tanh - = 1 (32) 

z z 

implying 

qy F ml.26n (33) 
Putting this value back in eq. (31) we obtain 

5 p 2 ^0.754 A BCS (34) 

From the expansion of the gap equation around 5p 2 it is easy to obtain 

A 2 « 1.76 5/i 2 (5p 2 -Sfi) (35) 

According to ref. [10] the difference between the grand potential in the 
superconducting state and in the normal state is given by 

^LOFF — ^normal = ~ [ (36) 

JO 9 

Using eq. (18) and eq.(34) we can write 

dg p dA BCS p d8p 2 



g 2 2 A BC s 2 5p 2 



(37) 



Noticing that A is zero for Sp 2 = 5p we are now able to perform the integral 
(36) obtaining 

^LOFF - ^normal ~ ~0.44 p{5p - 5 p 2 f (38) 
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We see that in the window between the intersection of the BCS curve and 
the LOFF curve in Fig. 2 and S/j,2 the LOFF phase is favored. Furthermore 
at the intersection there is a first order transition between the LOFF and 
the BCS phase. Notice that since 5ji2 is very close to 5fi± the intersection 
point is practically given by 5[ii. In Fig. 2 we show also the behaviour of the 
condensates. Altough the window (Sfxi,5fi 2 ) — (0.707, 0.754) A BC s is rather 
narrow, there are indications that considering the realistic case of QCD [19] 
the windows may open up. Also, for different structures than the single plane 
wave there is the possibility that the windows opens up [19]. 

3 Crystalline structures 

The ground state in the LOFF phase is a superposition of states with different 
occupation numbers (N even) 

\0)loff = £>|JV) (39) 

N 

Therefore the general structure of the condensate in the LOFF ground state 
will be 

mx)^(x)} = J2 c *N c N + 2e 2i " N - g {N\i;(x)iP(x)\N + 2) 

N 

= J2 A *e 2l « N - £ (40) 

N 

The case considered previously corresponds to all the Cooper pairs having the 
same total momentum 2q. A more general situation, although not the most 
general, is when the vectors <]n reduce to a set defining a regular crystalline 
structure. The corresponding coefficients A^ (linear combinations of subsets 
of the A^r's) do not depend on the vectors since all the vectors belong to 
the same orbit of the group. Furthermore all the vectors ^ have the same 
lenght [20] given by eq. (33). In this case 

(0|^(a;)^(a;)|0) = A^e 2 ^ (41) 

i 

This more general case has been considered in [12, 20] by evaluating the 
grand-potential of various crystalline structures through a Ginzburg-Landau 
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expansion, up to sixth order in the gap [20] 



Q = aA 2 + —A 4 + -A 6 (42) 
2 3 

These coefficients can be evaluated microscopically for each given crystalline 
structure. The procedure that the authors of ref. [20] have followed is to 
start from the gap equation represented graphically in Fig. 3. Then, they 



Figure 3: Gap equation in graphical form. The thick line is the exact prop- 
agator. The black dot the gap insertion. 

expand the exact propagator in a series of the gap insertions as given in 
Fig. 4. Inserting this expression back into the gap equation one gets the 



-•-o-#- 



Figure 4: The expansion of the propagator in graphical form. The darker 
boxes represent a A* insertion, the lighter ones a A insertion. 

expansion illustrated in Fig. 5. On the other hand the gap equation is 
obtained minimizing the grand-potential (42), i.e. 

a A + (3 A 3 + 7 A 5 + • • ■ = (43) 

Comparing this expression with the result of Fig. 5 one is able to derive the 
coefficients a, (3 and 7. 

In ref. [20] more than 20 crystalline structures have been considered, eval- 
uating for each of them the coefficients of Eq. (42). The result of this analysis 
is that the face-centered cube appears to be the favored structure among the 
ones considered (for more details see ref. [20]). However it should be noticed 
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Figure 5: The expansion of the gap equation in graphical form. Notations 
as in Fig. 4. 

that this result can be trusted only at T = 0. In fact one knows that in the 
(S/j,, T) phase space, the LOFF phase has a tricritical point and that around 
this point the favored crystalline phase corresponds to two antipodal waves 
(for reviews see [21, 22]). Therefore there could be various phase transitions 
going down in temperature as it happens in the two-dimensional case [23]. 

4 Phonons 

QCD at high density is conveniently studied through a hierarchy of effec- 
tive field theories. The starting point is the fundamental QCD lagrangian. 
The way of obtaining a lower energy effective lagrangian is to integrate out 
high-energy degrees of freedom. Polchinski [24] has shown that the physics 
is particularly simple for energies close to the Fermi energy. In particular he 
has shown that all the interactions are irrelevant except for a four-fermi inter- 
action coupling pair of fermions with opposite momenta. This is nothing but 
the interaction giving rise to the BCS condensation. This physics can be de- 
scribed using the High Density Effective Theory (HDET) [25]. In this theory 
the condensation effects are taken into account through the introduction of 
a Majorana mass term. The degrees of freedom are quasi-particles (dressed 
fermions), holes and gauge fields. This description is supposed to hold inside 
a shell around the Fermi surface of thickness 25 with the cutoff 5 smaller than 
the Fermi energy but bigger than the condensation energy A, A <C 5 <C Ep. 
Going at energies (measured with respect to the Fermi energy) much smaller 
than the gap energy A, all the gapped particles decouple and one is left with 
the low energy modes as Goldstone bosons, ungapped fermions and holes and 
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massless gauged fields according to the breaking scheme. The corresponding 
effective theory in the Goldstone sector can be easily formulated using stan- 
dard techniques. In the case of CFL and 2SC such effective lagrangians have 
been given in refs. [26] and [27]. The parameters of the effective lagrangian 
can be evaluated at each step of the hierarchy by matching the Greens func- 
tions with the ones evaluated at the upper level. The same technique will 
be applied in this Section to evaluate the propagation parameters of the 
Goldstone bosons (phonons) associated to the breaking of translational and 
rotational symmetries in the LOFF phase. The number and the features of 
the phonons depend on the particular crystalline structure. We will consider 
here the case of the single plane- wave [28] and of the face-centered cube [29] . 
We will introduce the phonons as it is usual for NG bosons [28], that is as 
the phases of the condensate. Considering the case of a single plane-wave we 
introduce a scalar field <&(x) through the replacement 

A(x ) = exp 2 ^* A -> e^ (x) A (44) 

We require that the scalar field <&(x) acquires the following expectation value 
in the ground state 

($(x)) = 2q-x (45) 

The phonon field is defined as 

-(f)(x) = $(ar) - 2q- x (46) 

Notice that the phonon field transforms nontrivially under rotations and 
translations. From this it follows that non derivative terms for (f>(x) are not 
allowed. One starts with the most general invariant lagrangian for the field 
<&(x) in the low-energy limit. This cuts the expansion of $ to the second 
order in the time derivative. However one may have an arbitrary number of 
space derivative, since from eq. (45) it follows that the space derivatives do 
not need to be small. Therefore 



J phonon 



$ 2 + ^c fe $(V 2 ) fc $ J (47) 



Using the definition (46) and keeping the space derivative up to the second 
order (we can make this assumption for the phonon field) we find 

Vonon = \ (0 2 - vlV ± (f) • V ± - Vj|0 • V||0) (48) 
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where 

— * 

V[|=n(n-V), V± = V-V[|, n = (49) 

We see that the propagation of the phonon in the crystalline medium is 
anisotropic. 

The same kind of considerations can be made in the case of the cube. The 
cube is defined by 8 vectors <f, pointing from the origin of the coordinates to 
the vertices of the cube parameterized as in Fig. 6. 




Figure 6: The figure shows the vertices and corresponding coordinates of the 
cube described in the text. Also shown are the symmetry axis. 

The condensate is given by [20] 

8 3 

A(x) = Aj2e mk '* = A Yl <? mHXi (50) 

k=l i=l,( £i =±) 

We introduce now three scalar fields such that 

(&®(x)) = 2\?\xi (51) 

through the substitution 

3 

A(x) — > A e ieMl){x) (52) 

t=l,(ei=±) 
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and the phonon fields 



= $®(x) -2\q\xi (53) 

Notice that the expression (52) is invariant under the symmetry group of the 
cube acting upon the scalar fields <&W(x). This group has three invariants 
for the vector representation 

I 2 {X) = \X | 2 , I 4 (X) = x\x\ + x\x\ + xjxl 

J 6 (X)=X 2 XfXf (54) 

Therefore the most general invariant lagrangian under rotations, translations 
and the symmetry group of the cube, at the lowest order in the time deriva- 
tive, is 

-^phonon ^ / , \ ^ ) 

i=l,2,3 

+ L S (/ 2 (V$W),/ 4 (V$«),/ 6 (V$«)) (55) 

Expanding this expression at the lowest order in the space derivatives of the 
phonon fields one finds [29] 

i=l,2,3 i=l,2,3 

- \ E W W ) 2 -c E ^ ( ^^ 0) (56) 

i=l,2,3 i<j=l,2,3 

The parameters appearing in the phonon lagrangian can be evaluated 
following the strategy outlined in [30, 31]. One starts from the QCD la- 
grangian and derives an effective lagrangian describing fermions close to the 
Fermi surface, that is at momenta such that p ~ pf but p — Pf ^ A. The 
relevant degrees of freedom are the fermions dressed by the interaction, the 
so-called quasi-particles [25]. Going closer to the Fermi surface the gapped 
quasi-particles decouple and one is left with the light modes as NG bosons, 
phonons and un-gapped fermions. It is possible to derive the parameters 
of the last description by the one in terms of quasi-particles evaluating the 
self-energy of the phonons (or the NG bosons) through one-loop diagrams 
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due to fermion pairs. The couplings of the phonons to the fermions are ob- 
tained noticing that the gap acts as a Majorana mass for the quasi-particles. 
Therefore the couplings originate from the substitutions (44) and (52). In 
this way one finds the following results: for the single plane-wave 



2 

and for the cube 




2 




(58) 



5 Astrophysical consequences 

A typical phenomenon of the pulsars are the glitches (for a review see [32]), 
that is sudden jumps in the period of the star. If pulsars are neutron stars 
with a dense metallic crust, the effect is explained assuming that some an- 
gular momentum is stored in the vortices present in the inner neutron su- 
perfluid. When the period of the star slows down , the vortices, which are 
pinned to the crystalline crust, do not participate in the slowing-down until 
they become unstable releasing suddenly the angular momentum. Since the 
density in the inner of a star is a function of the radius, it results that one 
has a sort of laboratory to study the phase diagram of QCD at zero temper- 
ature, at least in the corresponding range of densities. A possibility is that 
one has a CFL state as a core of the star, then a shell in the LOFF state and 
eventually the exterior part made up of neutrons. Since in the CFL state the 
baryonic number is broken there is superfluidity. Therefore the same mech- 
anism explained above might work with vortices in the CFL state pinned 
to the LOFF crystal. This could reinforce the ideas about the existence of 
strange stars (made of up, down and strange quarks). Considering the very 
narrow range of values for 5/j, in order to be in the LOFF phase one can ask 
if the previous possibility has some chance to be realized. Notice that using 
the typical LOFF value 5/i ~ .75Abcs one would need values of 5/i around 
15 4- 70 MeV. Let us consider a very crude model of three free quarks with 
M u = M d = 0, M s 7^ 0. Assuming at the core of the star a density around 
10 15 g/cm 3 , that is from 5 to 6 times the saturation nuclear density 0.15 x 10 15 
g/cm 3 , one finds for M s ranging between 200 and 300 MeV corresponding 
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values of 8/j, between 25 and 50 MeV. We see that these values are just in 
the right range for being within the LOFF window. 
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